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SKIN-STIFFENER COMBINATIONS 


By George E. Griffith and Georgene H. Miltohberger 


SUMMARY 


Temperatures and thermal stresses in typical skin-stiffener combi- 
nations of winglike structures subjected to aerodynamic heating have 
been obtained with the aid of an electronic differential analyzer. 
Variations were made in an aerodynamic heat-transfer parameter, in a joint 
conductivity parameter, and in the ratio of skin width to skin thickness. 
The results, which are presented in nondimens ional form, indicate that 
decreasing the joint conductivity parameter lowers both the interior 
and the average temperature ratios, increases the peak thermal stress 
ratios in the skin, and may considerably increase the peak stiffener 
stress ratios; increasing the aerodynamic heat-transfer parameter 
decreases the interior and average temperature ratios, increases the 
peak skin stress ratios somewhat, but greatly increases the peak stiff- 
ener stress ratios; and increasing the ratio of skin width to skin 
thickness produces only moderate decreases in the peak skin stress 
ratios while moderately increasing the peak stiffener stress ratios. 


INTRODUCTION 


Transient thermal stresses produced by aerodynamic heating of super- 
sonic aircraft may result in buckling, warping, flutter, loss of stiffness, 
or other deleterious effects which might seriously penalize the performance 
of the aircraft. In order for the designer to eliminate or minimize these 
effects, he must know how the thermal stresses are affected by such factors 
as flight conditions, materials, and the presence of structural disconti- 
nuities such as joints. At present, little information is available con- 
cerning the effects of joint conductivity on the thermal stresses; such 
information is needed for the alleviation of these thermal stresses and 
their accompanying undesirable effects. 

The purpose of the theoretical investigation described herein was to 
determine the effects of joint conductivity on the temperatures and the 
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thermal stresses in skin-stiffener combinations of winglike structures 
for different aerodynamic conditions. The work was carried out in non- 
dimensional form with the aid of an electronic differential analyzer. 
Geometrical effects were included by varying the ratio of skin width 
to skin thickness, aerodynamic-heating effects by varying the ratio 
of the product of the boundary-layer heat-transfer coefficient and the 
skin thickness to the coefficient of thermal conductivity of the material, 
and joint-conductance effects by varying the ratio of the product of the 
joint heat-transfer coefficient and the skin thickness to the coefficient 
of thermal conductivity of the material. 

In the present paper the variations in the parameters just described 
are discussed in terms of their influence on the maximum, mi ni mu m, and 
average temperatures, on the temperature drop across the joint, and on 
the maximum stresses in the skin and stiffener. These effects are illus- 
trated by an example and then discussed in detail. Equations are 
developed in appendixes for calculating the temperatures and stresses in 
the form of dimensionless ratios, and a discussion of the nondimensional 
parameters is also included as an appendix. 


SYMBOIS 


^n 


w 


E 

f 

h 

h j 

h jV k 

ht s /k 

hr/cwtg 


cross-sectional area of element n 
width of skin 
depth of stiffener 

specific heat of material 
modulus of elasticity 
length of joint 

boundary-layer heat-transfer coefficient 
joint heat -transfer coefficient 

joint heat-transfer parameter 

aerodynamic heat-transfer parameter (Biot number) 
time parameter 
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k 

M 

*8 

"^W 

T 

T av 

T m 

ATj 

T 0 


w 

a 

a 


thermal conductivity of material 

Mach number 
skin thickness 

stiffener thickness 
temperature 

average temperature of skin-stiffener structure 
adiabatic -wall temperature 
temperature drop across joint 
initial temperature 

specific weight of material 
coefficient of thermal expansion 
stress 


0 

Ea(T AW - T 0 ) 


stress ratio 


0 


time 

T - T 0 

temperature ratio. — 

T AR - T o 


Subscripts : 
s skin 


w 


stiffener 
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THEORETICAL INVESTIGATION 
Problems Investigated 


The present investigation consisted of finding the temperatures 
and thermal stresses in skin-stiffener combinations which were symmetri- 
cal about the stiffener center line (fig. 1 ) and which were symmetrically 
heated along both skins . Because of symmetry of both the structure and 
the heating, only one-half the geometry, as shown in figure 1 , need be 
considered in the analysis. This geometry was considered to be at a uni- 
form initial temperature and was then subjected to a step-function change 
in aerodynamic-heating conditions. 

The structure was assumed to be heated in such a way that the heat- 
transfer coefficient of the boundary layer and the adiabatic-wall tem- 
perature did not change with time and were constant along the skin of the 
structure considered. S imi larly, the thermal conductivity of the joint 
was assumed to have a constant value. The effects of radiation, both 
internal and external, and the temperature variation through the skin 
thickness were neglected. It was assumed that there were no heat losses 
in the system: the only external transfer of heat occurred between the 

boundary layer and the skin; internal heat transfer took place by con- 
duction; and the heat transfer across the joint was considered to be 
dependent upon the joint conductance. Material properties were assumed 
not to change with temperature. In calculating the thermal stresses, 
the slight dissymmetry about an axis perpendicular to the skin and passing 
through the midpoint of the joint was neglected, and all deformations of 
the structure were assumed to be elastic. 

As previously mentioned, the work described herein was carried out 
in terms of nondimensional parameters associated with the particular type 
of problem considered in this investigation. Hence, the problems con- 
sidered can be defined in terms of certain of these nondimensional ratios. 
Changes in the geometry of the structure were achieved by varying the 
ratio of the skin width to the skin thickness ^b s j t s = 20, 40 , or 60) 
while all other geometric ratios were held constant (b^/tg = 20, 
tw/ts = 0.50, and f/t s = 3*875; the ratio f/t s = 3*875 was arbitrarily 
obtained by letting t s = 0.250, = 0.125, and f = 2.5 (ty + t s ) + 1/32) . 

The boundary-layer aerodynamic heat-transfer coefficient was contained in 
an aerodynamic heat-transfer parameter htg/k, commonly referred to as the 

Biot number; values of ht s y/k used were 0 . 02 , 0 . 2 , and 2 . The heat- 
transfer coefficient of the joint appears in the parameter hjtg^k. Prob- 
lems were solved for arbitrarily chosen values of hjt^k = 0.03, 0.15, 

O.36, and 00. The special case of hjt g /k = 0 (zero joint conductivity) 
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was solved analytically, as shown in appendix A, without using the elec- 
tronic differential analyzer. 

Thus it can he seen that the problems investigated can be described 
in terms of a geometric ratio b s /t s , an aerodynamic heat-transfer param- 
eter htg/k, and a joint conductivity parameter hjt^/k. Problems were 

solved for all combinations of the parameters mentioned in the preceding 
paragraph, except that for b s /t s = 60 and ht^/k = 2 the only problems 

solved were for hjt s /k = 0 and °o. 


Theoretical Solution 

Inasmuch as the determination of the thermal stresses was the pri- 
mary objective of this investigation, and since their determination 
depends upon the distribution of the temperatures, a prerequisite to 
finding the stresses is a knowledge of the temperatures. Once the tem- 
peratures are known, the thermal stresses can be readily obtained (see, 
for example, ref. 1 ) ; however, the evaluation of the stresses, although 
straightforward, can be cumbersome - depending upon the problem, the 
simplifying assumptions, and the method of solution. 

Solutions to the well-known partial differential equation describing 
the conductive heat flow in a body are difficult to evaluate by ordinary 
means for the problem under consideration. Recourse is usually made to 
a mathematically simpler description of the problem, generally by dividing 
the structure into a number of elements of finite size (often referred to 
as lumps). An ordinary differential equation, continuous in time, can 
then be used to describe the thermal equilibrium of each element, the net 
result being a system of simultaneous ordinary differential equations 
equal in number to the number of elements. A closed-form solution for 
the temperature of any element at any time (provided there are several 
elements) can only be obtained through a cumbersome and tedious procedure . 
However, temperature histories for the individual elements can readily be 
obtained by using either a digital or an analog computer. In this inves- 
tigation solutions were obtained with the aid of an analog computer known 
as an electronic differential analyzer. In order to check the solutions 
for the temperatures a numerical procedure was used to solve sets of 
finite-difference equations for a few problems wherein the geometry was 
approximated by using fewer elements and the joint conductivity was con- 
sidered infinite. 

Temperatures . - For a given problem, the cross section of the struc- 
ture was divided into 15 elements, such as shown in figure 2. (Fig. 2 is 
a scale drawing of the divided geometry for b s /t s = 40.) The elements 
varied in size so that more elements could be concentrated in the region 
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of the joint and so that, at the same time, the number of elements could 
be kept to a reasonable minimum. An ordinary differential equation was 
used to describe the thermal equilibrium of each element, the net result 
being a system of 15 simultaneous ordinary differential equations. The 
development of these equations is given in appendix B, and the final set 
of these equations, written in nondimens ional form, is shown in matrix 
notation at the end of this- appendix. The following differential equation 
for element 5 is presented because it contains all the types of terms 
found in the complete set of equations . In nondimensional terms , the heat 
balance for element 5 is given by 


d fb - T o\ + - T o ) + 

- Tq/ P V t aw - T q/ 



C l^ fe - Tq V P- 1 foil - T 0 N 
p Vm - T o J _k_ + 1 + Jw_\T AW - t 0 

Vs 2 2t s 


where 

0 = hr/ cwt s 

P = ht s/ k 


and C]_i, C 12 , and are dimensionless geometric ratios similar to 

those shown for equation (B3) in appendix B; the subscripts used are the 
same as those in appendix B. The temperature ratio for any element is 
simply the ratio of the temperature rise of the element above the Initial 
temperature to the maximum possible temperature rise. 

Thermal stresses .- The assumptions previously mentioned, that 
material properties do not change with temperature and that the deforma- 
tions of the structure are elastic, greatly reduce the labor of computa- 
tion for finding the thermal stresses. If, in addition, there is no 
thrust or resultant bending moment acting on the cross section (as 
assumed in the present case) , the stress at a point then is directly 
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proportional to the difference in the average temperature of the structure 
and the temperature at that point. At any point i, the nondimensional 
stress ratio (which is treated more fully in appendix C) is given by 


Ea ( T Atf " T o) 


T i - T 0 , T av ~ T 0 


T AW - T o 


T AW - T o 


where the stress at a point is given as a fraction of a fictitious or ref- 
erence thermal stress associated with the maximu m possible temperature 
rise. Hence, the stress ratio can be either positive or negative (indi- 
cating tension or compression, respectively) but will always have an abso- 
lute value less than 1. 

Electronic differential analyzer .- The system of equations for the 

temperature ratios (appendix B) and the corresponding stress-ratio equa- 
tions were solved simultaneously by an electronic differential analyzer. 

A characteristic of this device is that variations in parameters exceeding 
a range of about 100:1 -frequently require reprograming, and, for this rea- 
son, cases for which ht s /k = 0.002 were not solved although these cases 

were included in the original program. Results are obtained from the 
analyzer by plotting the answers as functions of time or the answers may 
be read directly . from voltmeters at discrete times. The various features 
of this and other types of high-speed computers or analog methods are 
discussed in several publications, as, for example, in references 2 to 5. 


RESULTS AND DISCUSSION 
Data Obtained 


Three nondimensional temperature ratios ( T - t o)/( t aw - T o) and eight 
nondimensional stress ratios a] Ea (TAN - To) were plotted against the 
time parameter hT^/cwt s for each problem solved by the electronic differ- 
ential analyzer. Temperature ratios were obtained for the maximum temper- 
ature (element 1), the minimum temperature (element 15) , and the average 
temperature, and stress ratios were obtained for elements 1, k, 6, 7, 10, 
12, 13, and 15. The time scale for both temperature and stress ratios 
varied from problem to problem. After the analog computer had solved a 
problem, the stress ratios were examined to find the approximate value of 
the time parameter at which the maximum stiffener stress ratio had reached 
its peak value; the problem was then rerun so that voltmeter readings 
corresponding to all temperature ratios could, if possible, be obtained 
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shortly "before, at, and shortly after the peak stiffener stress ratio. 
These readings provided the only complete temperature distributions 
obtained and also aided in establishing values of the time parameter 
for peak stress ratios and maximum temperature -ratio drops across the 
joint. 

The overall accuracy of the results presented was affected by sev- 
eral sources of error, the magnitudes of which were difficult to deter- 
mine precisely. However, the maximum error in the temperature ratios 
was estimated to be less than ± 0 . 02 , and in the stress ratios to be less 
than ±0.05. The probable errors should be somewhat lower. 


Illustrative Example 

In order to illustrate the type of results obtained from the present 
investigation, the results for a single illustrative problem will be dis- 
cussed before examining the overall findings. The results for the illus- 
trative problem will be interpreted in terms of a specific material, 
geometry, and aerodynamic boundary condition. However, the general 
results of the succeeding sections will be discussed in terms of nondi- 
mensional parameters. 

Consider, for illustrative purposes, an aircraft wing made of 
type 302 stainless steel and flying at an altitude of 25,000 feet at a 
Mach number of 0.9- After reaching an equilibrium temperature of 33° F 
(corresponding to' the adiabatic -wall temperature), the aircraft is 
suddenly accelerated to a Mach number of 3.0 and then maintains this 
speed. (See fig. 3(a).) The approximate aerodynamic conditions in 
this second stage of flight correspond to the following: 


T 0 = 33° F 

M = 3-0 

Altitude = 25,000 ft 
h = 0.025 Btu/(sq ft) (sec) (°F) 

T aw = 66 2° F 

The portion of the structure to be analyzed consists of one-half of a 
skin and stiffener combination which is symmetrical about the stiffener 
center line and is similar to that shown in figure 1. Thus, with the aid 
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of figure 1, the geometry can be completely defined with the following 
quantities : 


t g = 0.250 in. 
b s = 10.000 in. 
by. = 5*000 in. 
iV = 0.125 in. 

f = O .969 in. 


The material properties for type 302 stainless steel were taken as 
(see ref. 6) : 

c = 0.12 Btu/(lb)(°F) 

w = 0.29 Ib/cu in. 

k = 0.0026 Btu/(ft)(sec)(°F) 

E = 28.0 X 10 6 lb/sq in. 
a, = 9*6 X 10"^ in.j (in.)(°F) 

The value of the joint conductivity used for this example is 

hj = 67.4 Btu/(sq ft)(hr)(°F) 

= O.OI 87 Btu/(sq ft)(sec)(°F) 

This value of hj is lower than expected in the usual stainless -steel 

joint (see, for example, ref. 7 ) and corresponds to one in which some 
type of insulation was used. 

The conditions outlined correspond to the problem for which 

ht s /k = 0.20 
hjt s /k = 0.15 
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The results of the analog computation for the illustrative problem 
are shown in figures 3 and 4. Figure 3 shows the change with time in the 
maximum skin temperature (element 1 ), the minimum temperature in the 
stiffener (element 15 ), the average temperature , and the corresponding 
stresses for elements 1 and 15 . The stresses shown are the maximum, in 
absolute value, for either the skin or the stiffener. Whereas , as seen 
in figure 3 (b), the skin temperature rises rapidly and approximates the 
maximum possible temperature T^y in about 4 minutes, the stiffener tem- 
perature for element 15 lags considerably and requires a considerably 
longer time to approach T^y. The average temperature of the cross sec- 
tion is much closer to the maximum skin temperature than to the minimum 
stiffener temperature. 

The difference between the average temperature and the skin or 
stiffener temperature is indicative of the magnitude of thermal stress 
experienced by element 1 or 15 . The difference between the average and 
the stiffener temperatures is considerably greater than the difference 
between the average and the skin temperatures ; this is reflected in the 
stress plots in figure 3(c). Since the temperature differences also have 
different signs, the stresses are of different character, the skin being 
in compression and the stiffener in tension, as expected. The peak values 
of the stresses occur at about 2 minutes and 3 minutes for the skin and 
stiffener, respectively. These times are only a small fraction of the 
time required for the entire structure to approach the equilibrium tem- 
perature but are fairly close to the time required for the skin tempera- 
ture to stabilize. 

Plotted in figure 4 are the^ temperature and stress distributions for 
three different times: 0.4 minute, 1.1 minutes, and 1.9 minutes. For 

simplicity the temperatures and stresses in the skin to the right of the 
joint have been omitted; the results, if plotted, would be very similar 
to those shown for the skin to the left of the joint (represented by the 
symbols and solid lines between points A and B). Up to 1. 9 minutes, the 
skin temperatures have risen sharply, whereas the temperatures of ele- 
ments 14 and 15 in the stiffener have changed very little. The tempera- 
ture drop across the joint has increased slightly from 0.4 minute to 
1.1 minutes, but there is a slight decrease from 1.1 minutes to 1.9 min- 
utes; this indicates that the time of maximum temperature drop across 
the joint may be close to 1.1 minutes. (An independent but approximate 
check of the results gave a time of 1.1 minutes for the maximum tempera- 
ture drop across the joint for this problem.) As can be seen from fig- 
ure 4 , the temperature drop across the joint has a substantial effect on 
the stiffener-flange stresses. However, there is no apparent connection 
between the magnitude of the temperature drop across the joint and the 
magnitudes of the peak values of the maximum stresses in the skin or 
stiffener (elements 1 and 15 , respectively), nor is there any connection 
with the times at which these peak stresses occur. 
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The results shown, which encompass several simplifying assumptions 
(particularly in regard to the stresses, since, for example, the propor- 
tional limit has probably been exceeded), indicate that, although the 
type 302 stainless steel may well withstand the effects of the tempera- 
ture level, thermal stresses may develop which are so severe that they 
dictate design changes. 


Temperature and Stress Ratios 

Three different sets of geometric ratios were used in the present 
investigation. These ratios were assigned the following values: 

b s /t s = 20, 40, 60 
s = 20 

^w/^s = 0*50 

f/t s = 3.875 


Only the ratio of skin width to skin thickness was varied while the 
ratios of stiffener depth, stiffener thickness, and length of contact 
flange of the stiffener to the skin thickness were held constant. For 
the purpose of discussion it is appropriate to consider first the 
results which were obtained for a single set of geometric ratios; the 
geometry selected is for b s /t s = 20. Results are discussed for this 
geometry for three values of the aerodynamic heat-transfer param- 
eter htg/k (0.02, 0.2, and 2) and for five values of the joint con- 
ductivity parameter hjt s y/k (0, 0.03, 0.15, O.36, and «,). The results 
shown for the case of hjt g ^k = 0 were not calculated by the electronic 

differential analyzer but were calculated as shown in appendix A. For a 
brief discussion of all the nondimens ional ratios or parameters see 
appendix D. 

ht s /k = 0.02.- Figure 5( a ) shows the variation in the maximum, mini- 
mum, and average temperature ratios (T - Tq) j (T^ - Tq) plotted against the 
time parameter hr j cwt s for different values of the joint conductivity 
parameter hjt s /k. The aerodynamic heat-transfer parameter was kept con- 
stant and was the lowest value used in this investigation (htg/k = 0.02). 
The temperature ratios for element 1 of the skin for cases other than 
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that for which hjt s /k = 0 fell slightly "below the curve for no con- 
duction along the skin ^hjt g ^k = 0) hut otherwise showed no particular 

trend and lay within the narrow limits of the too solid curves shown. 

This small lack of consistency probably reflects inherent inaccuracies 
in the analog computer. The average temperature ratio and particularly 
the temperature ratio for element 15 decrease in value as hjt g ^k 

decreases , especially as hjtg/k approaches zero. For this particular 
set of conditions, the temperature -ratio changes are fairly moderate for 
values of the joint conductivity parameter between 0.03 «>• however, 

far greater decreases in the average temperature ratio and the tempera- 
ture ratio of element 15 of the stiffener occur between joint-conductivity- 
parameter values of 0.03 and 0. For some materials a value lower than 
hjt s ^k = 0.03 can probably be obtained only by designing the joint to 

prevent the flow of heat across it. For example, for h.^t s jk = 0.03 
and a value of hj = 100 Btu/(sq. ft)(hr)(°F), the following skin 
thicknesses would be required for the materials listed: titanium 

and Inconel X, 0.030 inch; type 302 stainless steel, 0.034 inch; SAE 
1020 steel, 0.143 inch; FS1 magnesium, 0.194 inch; and 7075 -T aluminum 
alloy, 0.250 inch. For some of these materials, then, the more moderate 
changes which occur between joint-conductance ratios of <» and 0.03 are 
indicative of what may be expected for uninsulated joints. 

Shown in figure 5(b) are the stress ratios corresponding to the tem- 
perature ratios for elements 1 and 15, the maxi mu m values for the skin and 
for the stiffener. Appreciable increases in the peak values of the stress 
ratios for both, the skin and the stiffener occur as the joint conductivity 
parameter decreases to zero, but for values of hjt s / k between °° and 

0.03 the increases are more moderate. The positive and negative signs of 
the stress ratios indicate, as expected, that the cool stiffener is in 
tension while the hot skin is in compression. It can be noted that 
values of the peak stress ratios, for either the skin or stiffener, occur 
at values of the time parameter which increase with decreasing values of 
the joint conductivity parameter and, also, that peak values for the skin 
and stiffener for a particular set of conditions do not occur at the same 
value of the time parameter. 

ht s /k = O.2.- Figure 6(a) shows the effect on the temperature ratios 

(maximum, minimum, and average) of varying the joint conductivity param- 
eter for the same geometry (b s /t s = 20) and an aerodynamic heat-transfer 

parameter of 0.2. The results show that there is ah insignificant change 
produced in the maximum temperature ratio (element 1 in the skin) , a 
noticeable lowering of the average temperature ratio when the joint con- 
ductivity parameter decreases in value, and a similar but more pronounced 
change in the minimum temperature ratio (element 15 in the stiffener). 
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Figure 6(b) shows the corresponding effect on the stress ratios. 

The peak values of the maximum skin- stress ratio (element 1) increase 
as the joint conductivity parameter decreases; this is due to the fact 
that while the skin temperature ratio did not change appreciably, the 
average temperature ratio decreased somewhat as the joint conductivity 
parameter decreased, thereby causing a corresponding increase in the 
difference in the temperature ratios (upon which the size of the stress 
ratio depends). The peak values of the stiffener stress ratios (ele- 
ment 15) are much larger than for the skin and tend to decrease very 
slightly with decreasing values of the joint conductivity parameter. 

Thus, for this particular combination of configuration and aerodynamic 
heat-transfer parameter, there is very little change produced in the 
peak stiffener stress ratio by varying the value of the joint conduc- 
tivity parameter. However, the results indicate, as might be expected, 
that the case of zero joint conductivity does not provide an upper 
limit to the peak stiffener stress ratio. 

ht s y/k = 2.- Similar results for both the temperature and stress 
ratios for the same geometry (b s /t s = 20) and an aerodynamic heat-tr ans fer 
parameter of 2 are shown in figure 7. This value of ht s /k corresponds 
to rather severe heating rates. Figure 7(a) shows that increasing hjt s /k 

from 0 to 00 produces no discern i ble change in the temperature ratios for 
elements 1 and 15 and only a moderate increase in the average temperature 
ratio. Hence, for this high value of ht g /k the effect of the joint con- 
ductivity parameter is relatively unimportant, since, as expected, the 
skin temperature ratio approaches its limiting value before heat can be 
conducted into the interior, regardless of the conductivity of the joint. 
This results, as seen in figure 7(b), in only moderate change in either of 
the stress ratios. The peak skin stress ratio increases as the joint con- 
ductivity parameter decreases. However, it can be noted that the peak 
stiffener stress ratio again decreases with decreasing h^tg^k; this is 

expected since the average temperature ratio decreases while the stiffener 
temperature ratio stays close to zero, so that the differences in the tem- 
perature ratios also decrease. Thus, it can again be seen that, at the 
larger values of ht s /k, the case of hjt^/k = 0 does not constitute an 

upper limit to the maximum stiffener stress ratio. 

b s /t s = 20.- The results discussed so far are grouped in figure 8 so 

that, for this particular geometry (b s /t s = 20 ) , the influences of the 
aerodynamic heat-transfer parameter and the joint conductivity param- 
eter can be seen simultaneously. Curves are included for values of 
hjt s /k = 0, 0.03, and co, so that both the change occurring between 
hjt s /k = 0 and 0.03 and the change between h^tg^k = 0.03 and <» can 
be observed. When the range of hjt s /k from 0 to 0.03 is excluded, the 
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results indicate that the aerodynamic heat-transfer parameter may produce 
more change in both the temperature and. peak: stress ratios than is pro- 
duced. by changing the joint -conductivity ratio from 0.03 to «>. (This, 
however, includes an exceptionally wide range in htg/k.) In general, 

even greater changes occur as hjt^/k diminishes from 0.03 to 0. 

b g /t s = 40 and 60.- Plotted in figures 9 and 10 are composite results 
similar to those of figure 8 but for values of the geometric ratio 
h s /t s = 40 and 60, respectively. Figure 10 shows no temperature- or 
stress-ratio curve for hjt s /k =0.03 and ht s /k = 2, since, for this 
geometry and htg/k = 2, the only case solved by the electronic differen- 
tial analyzer was for hjt s ^k = co. However, it is apparent from an inspec- 
tion of figure 10 that, as previously mentioned (see the section for 
ht s /k =2), the effect of varying hjt s /k is small when ht s jk = 2, espe- 
cially so for h s Jt s = 60. The results presented in figures 9( a ) and 10(a) 
show that the temperature ratios behave as expected; there is little or no 
change in the temperature ratio- for element 1, but both the average tem- 
perature ratio and that for element 15 decrease with increasing ht s /k 
and with decreasing hjt s /k. Note that for figures 8(a), 9(a), and 

10(a) the curves for element 15 for ht s ^k = 2 and for any case where 
hjt s yk = 0 lie close to or along the zero axis. 

The peak values of the skin stress ratio for element 1 increase with 
increasing ht^k and with decreasing h^t^k. (See figs. 9(b) and 

10 (b).) For element 15 the peak stiffener stress ratios also increase 
with increasing htg^k. However, these stiffener stress ratios may either 

increase or decrease with decreasing h^tg^k; this is so because when 
htg/k is small the difference in the average temperature ratio and that 
for element 15 always increases as h j t s jk decreases, whereas as ht g ^/k 

gets larger the difference in these temperature ratios begins to decrease 
with decreasing h^tg^k. The results shown in figures 8(b), 90 3 ), and 

10(b) again illustrate that, although the case of zero joint conductivity 
constitutes an upper limit to the peak skin stress ratio, the upper limit 
to the peak stiffener stress ratio is not that for h^t s yk = 0. 

Effect of geometric ratio b s /t s . - Although the overall effects of 

changing ^ s /t s - the only geometric change effected - can be seen by 

studying figures 5 to 10, some of these results have been replotted in 
figures 11, 12, and 13 in order that they may be more readily studied. 

These latter figures show the temperature and stress ratios for 
htg/k = 0.2, for values of hjt s /k = 0, 0.03, and respectively, and 

for values of b s /t s = 20, 40, and 60. 



NACA TN 3699 


15 


Figure 11 shows results for the case of zero joint conductivity 
(hjt s /k = 0) which were obtained as outlined in appendix A (and not by 
using the electronic differential analyzer) . The only change in the tem- 
perature ratios occurs in the average temperature ratio, which increases 
as t> S/ /t s increases. Figure 11 also shows that the stiffener stress 

ratios increase and that the skin stress ratios decrease as i> s /t s 
increases, with the major changes occurring between b s /t s = 20 and 40. 

The values of the time parameter at which the peak values of these stress 
ratios occur are theoretically infinite, but, for practical purposes, can 
be seen to have been reached at values of hr/cwt s of about 6 and 4 for ■ 

the stiffener and skin, respectively. In general, the major result of 
increasing the ratio b^tg is effectively to put a greater proportion 

of the mass in the skin, thereby raising the average temperature ratio 
while having little or no effect on the magnitudes of the temperature 
ratios for elements 1 and 15. Thus, the difference between the skin tem- 
perature ratio and the average temperature ratio will decrease as' b g /t s 

is increased, and the corresponding thermal-stress ratio for element 1 will 
likewise decrease. On the other hand, the' difference in the stiffe ne r 
temperature ratio for element 15 and the average temperature ratio will 
increase with increasing bg^tg, so that the corresponding stress ratio 

for element 15 will also increase. 

Figure 12 shows results for ht s /k =0.2 and hjt s /k = 0.03. Again, 

as expected, there is little change in the temperature ratios for ele- 
ments 1 and 15, but a progressive increase occurs in the value of the 
average temperature ratio as ~b s /t s is increased, with more change 

occurring between b s /t g = 20 and 40 than between 40 and 60. The stress- 
ratio results are similar to those shown in figure 11 (b) for hjt g y/k = 0 

and, again, there is a slightly larger increase in the stiffener stress 
ratio and a correspondingly larger decrease in the skin stress ratio 
between b s /t s = 20 and 40 than between 40 and 60. The value of the 
time parameter for which peak values of the stiffener stress ratios 
occur decreases as b s /t s increases, whereas there is apparently little 

change in the time parameter for peak values of the skin stress ratios. 


The results shown in figure 13 for ht s /k = 0.2 and 
have essentially the same characteristics as those just discussed for 


hjt s /k = 


ht, 


3 /k = 


0.2 and 


hjt s y/k = 0.03 


and shown in figure 12. The collective 


results shown in figures 11, 12, and 13 illustrate that, for a point 
approximately 10 times the skin thickness or greater away from a he at 
sink, conduction toward the heat sink will make little or no difference 
in the temperature rise at that point. The effects of heat conduction 
near the heat sink will be more noticeable for low values of ht g /k 

and for high values of hj t^k . 
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The results shown in figures 11, 12, and 13 illustrate the effects of 
changing b g /t g for ht g /k = 0.2 .and are typical of the results in 

general. A more complete study of the effects of changing b g /t s can 

he made hy examining the results shown in figures 5 to 10. 


Peak Stress Ratios and Temperature -Ratio 
Drop Across Joint 

Variation of skin stress ratio (element 1) with hjt g yk .- The peak 

values of the stress ratios for element 1 of the skin are plotted against 
hjtg^k in figure 14(a) . This plot shows that there is a separate curve 

for each combination of ht s /k and b g /t g . (There is no curve for 
ht g /k = 2 and b g /t g = 60, since, of this group, the only problem solved 
was that for which hjt g y/k = °°.) The results show that the peak skin 
stress ratio decreases as b g /t g increases, which indicates that, as the 

proportion of skin area to total cross-sectional area increases, the 
average temperature ratio gets closer to the skin temperature ratio. The 
stress .ratios do not change much as hj t s/ k decreases, except that as 

■hjt g ^k approaches zero the peak skin stress ratio may increase appre- 
ciably, especially for low values of ht g ^/k. 

The values of the time parameter which correspond to the peak skin 
stress ratios are plotted in figure 14(b) against hjt g ^k. The results 

show that the value of the time parameter is essentially independent of 
the geometry and hence that a single curve can be drawn for each value 
of ht S y/k. The time parameter increases with increasing ht s /k and also 
increases with decreasing hjt g y/k, especially as hjt^k approaches zero. 

Variation of stiffener stress ratio (element 15) with hjt g /k.- Peak 
stiffener stress ratios are plotted against hjt g ^/k in figure 15(a) . 

Again, there is no curve for ht^k = 2 and b g /t s = 50. In some cases 
for which ht g ^k = 2, the stiffener stress ratio may not have reached its 

absolute peak value, since these stress ratios were leveling off as the 
computer was stopped. (See, for example, figs. 8 and 9*) For this reason 
the peak values shown in figure 15(a) are less reliable (only, it is 
believed, by a small amount) for these cases than for the others. As shown 
for the skin, there is a separate curve for each combination of ht g /k and 

b g /t g . The plotted results indicate that the peak stiffener stress ratio 
increases as ht_/k increases from 0.02 to 0.2, and again increases, but 
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not nearly so much, as htg/k increases from 0.2 to 2. The magnitudes 
of these changes depend upon hjt g y/k; below a value of h^t^k of about 
0.15, these changes diminish to zero (at hjt g ^k = 0); above 
hjt s /k = 0.15, the change with ht s /k is essentially constant. The 
results also show that as hjt g /k approaches zero the peak stress 
ratio may either increase appreciably at low values of ht^k or decrease 
slightly at high values of 


ht s /k. 


Values of the time parameter corresponding to the peak stiffener 
stress ratios are plotted against hjt g /k in figure 15(b). Since, as 

explained in the preceding paragraph, there was some difficulty in a 
precise determination of the peak stiffener stress ratio for some cases 
for which ht g ^k =2, it was impossible to determine the corresponding 

values of the time parameter, and hence no data are presented for these 
cases. Within the limits of accuracy of the analog computer, the results 
shown appear to be essentially independent of b g ^/t g . These limited 

results indicate that the time parameter increases appreciably as ht^k 

increases from 0.02 to 0.2. The results also reveal that the time parame- 
ter increases very slightly as hjt g /k decreases, except at the very 

‘ w ----- 


low values of h^tg/k where the increase is more pronounced 


Variation of temperature -ratio drop across joint with h.t /k.- The 

effects of varying the joint conductivity parameter hjt g /k on the tem- 
perature and stress ratios for elements 1 and 15 have been discussed in 
the preceding' sections . Also of some interest is the effect on the 
temperature -ratio drop across the joint. Shown in figure 16, plotted 
against hjt g ^k, are the maximum values of the temperature -ratio drop 

across the joint (that is, from element 6 to element 12) and the values 
of the time parameter at which these maximum drops occur. The results 
shown for the temperature -ratio drop and the time parameter are similar. 
Varying the geometric ratio b g y/t g from 20 to 60 seems to have little 

or no effect, as expected; hence, a single curve can be drawn through the 
points for each value of ht^/k. Both the maximum temperature -ratio drop 

and the corresponding value of the time parameter are affected by ht g /k 
and hjt^k; both decrease in value as ht g /k decreases and as hjt g ^k 
increases. This result is expected, since either decreasing ht s /k or 
increasing hj t s J k tends to permit freer conduction of heat in the 
structure . 


Neither the maximum temperature -ratio drops across the joint nor the 
corresponding values of the time parameter have been correlated with the 
peak skin or stiffener stress ratios and their associated time parameters. 
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It is Relieved that so many factors affect the temperature- and stress - 
ratio distributions and the magnitudes of the peak stresses that no such 
correlation can readily he made. Essentially, the same peak stresses may 
he produced in many ways. For illustrative purposes, figure 17 shows 
three temperature- and stress-ratio distributions. (These results were 
taken from three different problems at values of the time parameter 
slightly less than that required for peak stiffener stress ratios.) 
Although rather insignificant changes have occurred in the peak stiffener 
stress ratio (point D of element 15) and not much change in the peak skin 
stress ratio (point A of element 1 ),- the joint conductivity parameter can 
be seen to have had a tremendous effect on both the temperature and stress 
ratios in the vicinity of the joint. 


CONCLUDING REMARKS 


A theoretical investigation was conducted to determine the influence 
of joint conductivity on the temperatures and thermal stresses in aerody- 
namically heated skin-stiffener combinations of winglike structures. 
Varying the ratio of skin width to skin thickness, the aerodynamic heat- 
transfer parameter, and the joint conductivity parameter indicated the 
following results . 

Increasing the ratio of skin width to skin thickness raised the 
average temperature ratio without notably changing the other tempera- 
ture ratios, moderately decreased the peak skin stress ratios while 
moderately increasing the peak stiffener stress ratios, and had no 
apparent influence on the values of the time parameters associated with 
the peak stress ratios. 

Increasing the aerodynamic heat-transfer parameter produced the 
following results: considerably altered the distribution of the tempera- 

ture ratios by increasing the temperature -ratio drop across the joint, 
lowering. the minimum temperature ratio in the stiffener, and also lowering 
but to a considerably lesser extent, the average temperature ratio; appre- 
ciably increased the peak stiffener stress ratios while only moderately 
increasing the peak skin stress ratios ; and increased the values of the 
time parameter for both the peak stress ratios . 

Decreasing the joint conductivity parameter produced the following 
results: increased the termperature-ratio drop across the joint and 

lowered the average temperature ratio and the minimum stiffener tempera- 
ture ratio; increased the peak skin stress ratios appreciably only at low 
values of both the joint conductivity parameter and the Biot number; 
increased the peak stiffener stress ratio appreciably at the low values of 
the heat-transfer parameter but decreased the peak stiffener stress ratios 
slightly at the high values of the heat-transfer parameter; and had little 
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effect on the time parameter associated with the peak stresses, except 
to increase the time parameter sharply as the joint conductivity param- 
eter approached zero. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 14, 1956. 
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APPENDIX A 

SPECIAL CASE OF ZERO JOINT CONDUCTIVITY 
Temperatures 


For the case of zero joint conductivity, the temperature rise any- 
where in the skin is given hy the well-known differential equation 
(normally used when conduction effects are ignored) 

pm 

CWt s IF = h(T AW “ T s) (A 1 ) 


The solution to equation (Al) can he mitten in terms of the nondimensional 
temperature rise of the skin (provided that, as previously assumed, c, w, 
and h are constant with time) 


T„ - T, 


0 _ 


t AN " T 0 


= 1 - e 


hT 

CWte 


(A2) 


Since no heat can he conducted across the joint, the entire' stiffener 
remains at the initial temperature; hence, the nondimensional temperature 
ratio for the stiffener is zero. Thus, the average nondimensional tem- 
perature of the structure is given hy 


^av T 0 
T M “ T 0 


15 

XI An 

n=l 




(A3) 


where A n is the cross-sectional area of element n. 


Stresses 

The thennal stress at any point in a structure depends upon the 
difference in the temperature at the point and the average temperature 
(see, for e xam ple, refs. 1 and 8); hence, for the skin the nondimensional 
stress ratio is 
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Ea(T AW - T 0 ) 


T s ~ + ^av ~ T 0 

T AW “ T 0 T AW " T o 





(Ah) 


and for the stiffener 


Ea ( T AW “ T 0) 


Vb 


15 


ta 


n=l 



(A 5 ) 
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APPENDIX B 

DIFFERENTIAL EQUATIONS FOR CALCULATING TEMPERATURES 


Equations of heat balance are derived for three typical elements in 
a s kin and stiffener combination such as shown in figure 2. The assump- 
tions made in this analysis regarding the temperatures are as follows: 

(1) The heat-transfer coefficient of the boundary layer is constant 
along the skin at any time. 

(2) The temperature is uniform throughout any element but for ana- 
lytical purposes is considered concentrated at the center of the element; 
end elements 1, 10, and 15 are considered to be one-half of symmetrical 
elements . 

(3) Each element receives heat directly only from adjacent elements 
and, if a skin element, from the boundary layer. 

(h) Radiation effects can be neglected. 

(5) Material properties do not change with temperature. 

(6) Joint conductivity does not change with temperature. 

( 7 ) The skin and stiffener are made of the same material. 


Typical Skin Element 

Consider first a typical skin element (other than element 5 or 6), 
such as element 3, with the quantities of heat which influence its tem- 
perature (per unit depth, where the depth direction is perpendicular to 
the cross section of the structure) shown by the following sketch: 


_ kt„ 

<*2 = 7-M?2 “ *3) 

L2,3 


QbL = ^(^AW " T 3^ 
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where 

Q, rate of heat flow or heat stored 

L, distance between centers of elements 2 and 3 

2,3 

distance between centers of elements 3 and 4 

Thermal equilibrium requires that the rate of heat stored Q,^ be equal 
to the Sim of the rates of heat entering and leaving, Q - RT| + Qg + Q4, so 
that 


dT* kt kt 

a,t sb aT ' I'Ljfl'M - T 3 ) + j^(T 2 - Tj) H: " ? 3 ) (M.) 

3,4 


Subtracting Tq from each temperature, dividing through by hL^(T^j.j - To), 
and converting to nondimens ional time [where 9 = fa T ) gives equation (Bl) 

V CW V 


in nondimensional form 


& ( t 3 ~ T 0 


kt s ~ T 0 


de\TAw - T 0 J hL 3 L 2 ^\T AW - T 0> 


+ 1 + 


^3^2,3 


— t s V T 3 - Tq \ _ kt s fa - T 0 \ = x 
hL 3 L 3 , 4y v 1 ' AW - T 0 j hL 3^,4\ P AW " T 0/ 


which can be simplified to 


where 




d ft 3 ~ T o 


d0 V t AW " T 0 . 
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and G^, Cg, and Cj are dimensionless geometric ratios which, are the 

same as those appearing at the end of this appendix where the equations 
for all the elements are given in matrix form. 


Note that although the h used herein pertaining to the skin is 
described as the boundary-layer heat-transfer coefficient, this coeffi- 
cient could, in substance, be an effective heat-transfer coefficient, a 
combination of the heat-transfer coefficient of the boundary layer and 
the thermal conductivity of some surface insulator divided by its thick- 
ness t. This effective heat-transfer coefficient, which would then 
replace h in the calculations , would be related to the latter by the 
equation 


r 

^effective 



© 


insulator 


Using ^effective wou 13- give accurate results, provided little heat is 

stored in the insulator. Thus, the results reported may be interpreted 
for either insulated or uninsulated skins. 


Element Affected by Joint 


Before deriving similar equations for the elements affected by the 
joint, it is desirable to develop an expression for the heat transfer 
between elements separated by the joint such that this expression con- 
tains both the heat-transfer coefficient of the joint and the thermal 
conductivity of the material. Consider, for example, elements 5 an( ^ H 
which are separated by a joint of negligible thickness but of finite con- 
ductivity such that a drop in temperature takes place across the joint, 
as shown in the following sketch: 


Joint 



T, °F 


T11 



Across the joint, since there is no storage of heat. 


(B4) 


^(T 5 - T a ) = hj(T A - T b ) = -Mt b - T X1 ) 
a 5 A ll 



N 
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Solving equations (Bk) for T A and T-g in terms of and Til gives 


h. 


T A = 


_ V^3 a 3 a ii 


J T 5 + An X H 


T n 


11 


hi hi , 

J + -J-+ k 


A5 A 11 A 5 A 11 


(B5) 


and 


T B - 


*5 


T 5 + 


V A 11 


a 5 a h. 


IT- 


11 


h 


1 hi 

J + T^~ + 


A5 A n *5*11 


(b6) 


But, since there is no loss or storage of heat. 


hj(T A - T B ) . h>(T 5 - III) 


(B7) 


where h' is a coefficient which contains both the thermal conductivity 
of the material, divided by appropriate distances, and the heat-transfer 
coefficient of the joint. Solving equation (B7) for h' by substituting 
T a and Tg from equations (B5) and (b 6), respectively, gives 


h' 



(B8) 


Equation (b 8) is the expected result and is analogous to the well-known 
thermal series -resistance equation. (See, for example, ref. 9-) 

Since A^ = t s /2 and Aqq = \jf 2 , equation (b 8) can be written 
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The value of h' is largely governed by the value of hj, so that, as 
might be expected, the terms tg/k and t^k could be omitted with very 
little loss in accuracy when, hj is small (that is, when its reciprocal 
is large and; h' is therefore small). However, the values of the 
k terms are more significant and cannot be neglected as hj increases 

in value (and its reciprocal decreases). 


How consider element 5 with the quantities of heat which affect its 
temperature rise, as shown in the following 1 sketch: , 


% 



kt s 

* 5,6 


(t 6 - T 5 ) 


< 


where the notation is similar to that used in obtaining equation (Bl) . 
Hence , 


dTc kt s 

cwt s L 5 = hL 5 (T M - T 5 ) + — MT4 - T 5 ) + 

dT L 4,5 


kt. 


l 5,6 


(T 6 - T 5 ) + h’L 5 (T 11 - T 5 ) 


(BIO) 


Proceeding as before, equation (BIO) can be written in nondimensional 
form as / \ 


h * -v + 1 + 


'12 


,-l 


P 


* + 1 + ^ 




hjt s 2 2t 


4 


C 13 


06 " 


,-1 


k + l + 


0ii = 1 


(Bll) 


hjt s 2 2t g 
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The dimensionless ratios C-q, ^12l } an ^- are the same as those used 

in the text and in the equations at the end of this appendix. 


Typical Stiffener Element 


For a typical element unaffected by heat from the boundary layer, 
for example element 14, the quantities of heat which affect its tempera- 
ture rise are shown in the following sketch: 



kt w 

L l3,l k 


( t 13 “ T Hf) 


QX4 = cwlyL-^ 



As before. 


cwt w L l4 


^14 

dT 


kt w 

L l3A k 


(T15 - T &) 


+ 


kk w 

l i4,15 


( t 15 - T iU 


(B 12 ) 


Equation (B12) can be simplified to give 



"21 


0 i5 + 


"22 


i 


14 


C23 


0i 5 = ° 


(B13) 


Again, the coefficients C 2 j _ , C 22 > and C 23 are the same as those given 
in the matrix of the complete set of equations. Similarly, equations 
can be written for the remaining elements. 
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Nondimens ional Equations for Aerodynamically 
Heated 15-Element Structure 

The result of writing the equations for the remaining elements can 
"be conveniently expressed in matrix notation as follows: 

[%mj j^mj “ | p n| 


a i,i = a io,io = ^ + 1 + y 

-^ = A + + n 

■ LL , J - L as p 

C 1 

%2 - a 10,9 = " T 

a. - C l4 

1l,12 0- 

°2 

1,1 " *9,10 = f 

C 15 

a 12,01 = — 

1,2 = 1,9 = Jr + 1 •'x 

a X2 12 = V" + + n 

°k 

*2,3 = *9,8 - p 

a, - °17 

^2,13 ” p 

C 5 

a 3,2 = a 8,9 “ p" 

c l8 

*13,12 = — 

a 3,3 * a 8,8 = ^ + 1 + y 

a _ a C 19 

13,13 ~ aif + ~p” 

c 7 

- a 8,7 - p 

®20 

a i3,l4 = X 

Co 

%3 = “7,8 = p" 

*14,15 = 

a M = a 7,7=i- +1 + 7T 

14 , 14 = i- + x 

a 4,5 = 1,6 - ^ 

_ °23 

*14,1 = X 

1,4 = “6,7 = X 

C 2 I 4 . 

a 1,l4 = “ X 

“5,5 = a 6,6 = S- +1 + X + n 

a _ a C 2l 

a !5,i5 - de + x 

c 13 

1,6 = “6,5 = ~f 

All other a^ m = 0 


1,11 - ®6,i2 ~ a ii,5 = a i 2,6 = -n 
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where 


P 1 - p 2 - * * * " p io " 1 

P, = p = . . . = p =0 
11 12 15 

3 = ht s ^k 
0 = hr^ewtg 

0- ^ 

k + 1 + 
hjtg + 2 + 2tg 


and 

C 1 ~ t s 2 /hJsi } 2 

C 9 — c 8 “ c io 

C 17 = " t s 2 /L 1 2 L 12,13 

c 2 = ~ t s 2 f L 2^L J 2 

c 10 = _t s 2 /L4 L 4,5 

c l8 = - t s 2 /l a .3 L i2,13 

— t ■ 

0 

CVJ 

0 

1 

11 

tn 

0 

C 11 - -tsf 151^,5 

c 19 = “ c l8 “ c 20 

C 4 t s^/ L 2L 2 ^ 

C 12 = -c ll “ c 13 

C 20 = 

c 5 = _t s : 2 / L 3L2,3 

C 13 = - t s 2 / I 5 I 5,6 

C 21 = “ t s 2 /l'i4 L 13,l4 

c 6 = ~ c 5- “ c 7 

C l4 = t s 2 /LilLil,12 

C 22 = " c 21 “ c 23 

c 7 = -t s 2 / L 3 L 3,i}. 

C 15 = “ t s 2 / L l2 L ll J) 12 

c 23 = “ t s 2 / I G.4 L llt-,15 

C 8 = _t s 2 / L 4 L 3^ 

C l6 = " C 15 " C 17 

C 24 = ts2 / L i5 L L4,i5 
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APPENDIX C 

/ 

l 

THERMAL-STRESS EQUATIONS 


Since the stresses depend upon the temperatures, the assumptions 
used in calculating the temperatures also apply to the stresses . In 
addition, the following assumptions were incorporated in calculating 
the stresses: I 

(1) The cross section and temperature distribution are constant 
in the direction perpendicular to the cross section, the cross section 
under consideration is far enough from the end, of the structure to be 
free of any end effects, and only self -equilibrating thermal stresses 
act on the cross section. 

(2) The section considered is one-half of a .section which is symmet- 
rical about the stiffener center line. 

I 

(3) The slight dissymmetry about an axis .through the joint center 
line is considered to have a negligible effect and is thus ignored. 

(4) The stresses and strains are entirely elastic, and the modulus 
of elasticity does not change with temperature.. 

In accord with the aforementioned assumptions and those used in the 
determination of the temperatures , the stresses may be calculated directly 
(see refs. 1 and 8) from the equation 


°i — Ea(-T^ + T av ) 


= -EaT ± + 15 


15 - 
Ea y 
. n=l 


n=l 


*n 


(Cl) 


where 

rr . 

stress in element i 

1 

T i 

temperature of element 

T 

av 

average temperature of 


15 / 15 

H VWH An 

n=l / n=l 
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Equation (Cl) can be written in terms of a nondimensional stress 
ratio and nondimensional temperature ratios as 


^(TflW “ T 0 ) 


-13 rp _ ro 

is — ±Q_ * 

T i - T o t aw - T o 

t A w “ t o 1 5 

XI A n 

n=l 


15 




n=l 




15 


n=l 




(C2) 
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APPENDIX D 

NONDIMENSIONAL RATIOS 


The problems solved in this investigation and the results obtained 
have been described in terms of the natural parameters or nondimensional 
ratios associated with the particular type of problem investigated. (See 
appendixes B and C . ) The parameters which define a particular problem, 
once they have been assigned specific values, are t> g /t g , t^/tg, f/t g , 

t^/tg, htgyte, and hj-tg^k. Of these parameters, the first four define 

the geometry and the last two the boundary conditions. For prescribed 
values of these six nondimensional ratios (together with initial condi- 
tions) there exists a unique set of nondimensional temperature 
ratios (T - T 0 ) /(%.j - T 0 ) and stress ratios cr/Ea(T AW - T 0 ) which vary 

with the time parameter hr j cwt g (also nondimensional). 


Physical Significance 

The parameters b g ^t g , h^tg, f/t g , and t^tg are obvious geo- 
metric ratios which give the dimensions of the structure in terms of the 
skin thickness. The parameter ht s ^k, the Biot number, is the ratio of 

the surface heat-transfer coefficient to the internal heat-transfer coef- 
ficient (the coefficient of thermal conductivity of the material), multi- 
plied by a basic dimension which is characteristic of this problem, the 
skin thickness. Similarly, the parameter hjt g /k is the ratio of the 

joint heat-transfer coefficient to the internal heat-transfer coefficient, 
multiplied by the skin thickness. The ratio (T - T 0 )/ (T^W " ^0) < 3 - e fi nes "the 

temperature rise (above the initial temperature ) at a point as a fraction 
of the maximum possible temperature rise. The stress ratio 0/Ea(T^y - Tq) 

gives the stress at a point in terms of a reference thermal stress asso- 
ciated with the maximum possible temperature rise (T^j - Tq). The time 

parameter hr/cwt s is the ratio of the surface heat-transfer coefficient 

to the heat absorptive capacity of the skin, multiplied by the time, or' 
the time parameter may also be considered as the ratio of the actual time 
to the time constant of the skin (See eq. (A 2 ) in appendix A). 


Application to Specific Conditions 

For a given problem (a prescribed geometry or a set of geometric 
ratios and prescribed boundary and initial conditions), the results (the 
variation of the temperature and stress ratios with the time parameter) 
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may be applied to any material for various combinations of geometric size, 
aerodynamic heat-transfer coefficient, and joint conductivity. Once the 
type of material has been designated, the values of k, c, w, E, an d a 
become fixed. Then, selection of a value of h determines the value of 
t g or selection of t g determines the value of h, as dictated by the 

aerodynamic heat-transfer parameter ht s /k. The selected value of h can 
usually be attained under many aerodynamic conditions; these aerodynamic 
conditions must then be used to find both T AW and T 0 ; the actual values 

of the temperatures may then be converted from the nondimensional tempera- 
ture ratios. The values already selected for t g and k dictate the 

value of the joint heat-transfer coefficient to be obtained from the 
ratio hjt^k. Men combined with the set of geometric ratios, the value 

of t g determines the actual size of the structure. The values of T AW 
and Tq, which are dependent upon the aerodynamic conditions, and the 

material properties establish the magnitudes of the stresses (converted 
from the stress ratios). The times for the converted values of the tem- 
peratures and stresses depend upon the aerodynamic heat-transfer coeffi- 
cient, the material properties (c and w), and the skin thickness and 
are found from the time parameter hT/cwt g . Thus, the nondimensional 

results for any one prescribed problem may be converted to many different 
results depending upon the selection of the aerodynamic heat-transfer 
coefficient, the material, the skin thickness, the joint conductivity, or 
the time. Any one or a combination of these quantities may dictate the 
magnitudes of the temperatures and stresses and the boundary conditions 
for an actual structure. 
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Figure 1.- Dimensions and symbols for one-half of symmetrical 

skin- stiffener geometry . 



Figure 2.- Division of skin-stiffener geometry for analog calculations 
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Figure 3 *- Flight history and temperature and stress histories 

for illustrative problem. 
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(b) Stresses, 


Figure 5.- Concluded 
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(a) Temperatures. 


Figure 6.- The- effect on the temperatures and stresses of 


joint conductivity parameter 


hjt s /k. 


j/ts = 


20 and 


varying the 
ht s /k =0.2. 
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(b) Stresses. 


Figure 6.- Concluded. 
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(a) Temperatures. 

Figure 7 *- The effect on the temperatures and stresses of varying the 
joint conductivity parameter hjt g /k. b s /t g = 20 and ht g /k = 2. 
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(b) Stresses. 


Figure 7 .- Concluded. 
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(a) Temperatures. 

Figure 8.- The effect on the temperatures and stresses of varying 
hjt s ^k and ht s /k. b s /t s = 20. 








Pi f 





48 


NACA TN 3699 



O 2 4 6 8 10 


hr 

cwt s 


(a) Temperatures. 


Figure 10.- The effect on the temperatures and stresses pf varying 
hjt s ^k and ht s /k. h> s /t s = 60. 
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(b) Stresses. 


Figure 10.- Concluded 
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(a) Temperatures. 


Figure 11.- The effect on the temperatures and stresses of varying the 
geometric ratio "u s /t s . ht s /k = 0.2 and hjt s /k = 0. 
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(a) Temperatures. 


Figure 12.- The effect on the temperatures and stresses of varying the 
geometric ratio b s ^t s . ht s ^k =0.2 and hjtg^k = O.03. 
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(a) Temperatures. 

Figure 13 The effect on the temperatures and stresses of varying the 
geometric ratio b g y/t s . ht s ^k = 0.2 and hjt s /k = 00. 
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Figure 13.- Concluded 
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(a) Peak skin stress ratio. 


Figure 14. - Variation of peak skin stress ratio and time parameter with 

hjts/k. 
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(b) Time parameter. 


Figure Ik.- Concluded. 
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(a) Peak stiffener stress ratio. 


Figure 15.- Variation of peak stiffener stress ratio and time parameter 

with. hjt s ^k. 










